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The Double Chooz experiment has determined the value of the neutrino oscillation parameter ^13 
from an analysis of inverse beta decay interactions with neutron capture on hydrogen. This analysis 
uses a three times larger fiducial volume than the standard Double Chooz assessment, which is 
restricted to a region doped with gadolinium (Gd), yielding an exposure of 113.1 GW-ton-years. 
The data sample used in this analysis is distinct from that of the Gd analysis, and the systematic 
uncertainties are also largely independent, with some exceptions, such as the reactor neutrino fiux 
prediction. A combined rate- and energy-dependent fit finds sin^26'i3 = 0.097 ± 0.034 (stat.) ± 
0.034 (syst.), excluding the no-oscillation hypothesis at 2.0 a. This result is consistent with previous 
measurements of sin^2&i3. 

PACS numbers: 14.60.Pq 



Neutrino oscillations are well established in the three 
flavor paradigm and can be described by three mix- 
ing angles (0i2, 6*23, 6*13), a CP-violating phase 5, and 
two mass-squared differences {/S.m\i, Am§2)- Among the 
three mixing angles, ^13 is the smallest and has recently 
been revealed to be non-zero [lH3- The value of ^13 
is a critical input for plans to measure 5 and the neu- 
trino mass hierarchy. Furthermore, it may provide im- 
portant clues for physics beyond the Standard Model. 
The current best measurements of ^13 come from the re- 
actor j/g-disappearance experiments Double Chooz, Daya 
Bay, and RENO 0-0 • three experiments rely on 
the detection of the inverse beta decay (IBD) interac- 
tion, z/g -|- p — > -|- n, in Gd-doped liquid scintillator 
(LS). Typically these experiments search for a prompt 
positron signal followed by an 8 MeV gamma cascade 
from neutron capture on Gd. Background due to natu- 
ral radioactivity, which is predominantly below 4 McV, 
is largely suppressed. However, in Double Chooz it is 
also possible to search for a prompt positron followed by 
a 2.2 MeV gamma ray from neutron capture on hydro- 
gen, thanks to the low background environment in the 
detector. 

Though the latter analysis presents several challenges, 
it provides important benefits: a cross-check on the stan- 
dard Gd analysis and improved i>e energy spectrum shape 
information which is essential to our knowledge of ^13. 

In this letter we present an analysis of IBD interac- 
tions with neutron capture on hydrogen in the Double 
Chooz far detector. Following the same approach as in 
previous reports Q , this analysis compares the can- 
didate event rate and prompt energy spectrum shape to 
the Monte Carlo (MC) prediction. This analysis, how- 
ever, differs from those reported 0,0] in two major ways. 
First, the definition of the delayed signal is changed from 
the ^ 8 McV gamma cascade characteristic of a neutron 
capture on Gd to the 2.2 MeV gamma ray characteristic 
of a neutron capture on hydrogen. This change allows 



us to select a data set that is statistically independent of 
the Gd-based data set and has different systematic uncer- 
tainties and background characteristics. Second, because 
hydrogen captures occur in the undoped LS in addition 
to the Gd-doped region, a three times larger fiducial vol- 
ume is available for analysis. 

The Double Chooz far detector is located at a distance 
of ~ 1050 m from the two 4.25 GWth reactor cores of 
the Chooz Nuclear Power Plant, with a rock overburden 
of 300 meters water equivalent. The central region of 
the detector consists of three concentric cylinders, col- 
lectively called the inner detector (ID). The innermost 
cylinder is the 10.3 m^ target. This is surrounded by a 
7-catchcr (22.5 m-^). The target liquid is a PXE-based LS 
doped with Gd at a concentration of 1 g/1 while the 7- 
catcher liquid is an undoped LS. Outside the 7-catcher is 
the buffer, a 105 cm thick layer of non-scintillating min- 
eral oil contained in a stainless steel tank. Light from 
the target and 7-catcher volumes is collected by 390 low- 
background 10-inch PMTs installed on the inner wall of 
the buffer tank [§-ll|. Outside the buffer tank, and op- 
tically isolated from it, is the inner veto (IV), a 50 cm 
thick layer of liquid scintillator in a steel tank. The IV 
is equipped with 78 8-inch PMTs and serves as a veto 
for cosmic rays and fast neutrons entering the detector. 
The IV is surrounded by a 15 cm thick layer of demag- 
netized steel which suppresses 7-rays from radioactivity 
in the surrounding rock. Above the IV is the outer veto 
(OV) detector, a scintillator-strip-based muon tracking 
system. The OV system was installed during the data 
taking period, and about 2/3 of the data in this analysis 
benefit from OV use. A more detailed description of the 
entire detector can be found in Ref. 

The number of protons is estimated to be (6.747 ± 
0.020) X 10^9 in the target @ and (1.582 ± 0.016) x 
10'^'^ in the 7-catcher volume, the latter being based on a 
geometrical survey and measurements of the scintillator 
hydrogen fraction. 



3 



The IBD signal is a twofold coincidence of a prompt 
positron energy deposition, Eprompt, and a delayed 
gamma energy deposition, Edeiay, resulting from a neu- 
tron capture on hydrogen or Gd. The separation in time 
and space. At and Ar, of the coincident events are de- 
termined by neutron capture physics. Neutron capture 
times are 200 ^s in the 7-catcher and 30 fis in the target, 
where the presence of Gd greatly increases the neutron 
capture probability. In this analysis, where we search for 
Edoiay ~ 2.2 MeV without any fiducial volume cuts, we 
expect to detect candidates in both the target and 7- 
catcher. Given that only 13 % of the IBD interactions in 
the target volume are followed by neutron capture on hy- 
drogen [61] , 95 % of the signal events used in this analysis 
are located in the 7-catcher. 

Vertex reconstruction is based on a likelihood maxi- 
mization of the charge and timing of the pulses detected 
at each PMT It allows the spatial correlation of 
prompt and delayed events, effectively removing acciden- 
tal backgrounds. 

We reconstruct the energy of all events via two steps: 
(1) a total charge (Qtot) to photoelectron (PEtot) conver- 
sion; and (2) a PEtot to visible energy (Evis) conversion as 
done in the Gd analysis @ . The first step takes into ac- 
count a channcl-by-channel, non-linear gain calibration. 
The second step uses a light yield of ~ 230 PE/MeV, de- 
fined by the neutron capture peak on hydrogen in ^^^Cf 
calibration source data. By applying correction factors 
derived from spallation neutron data, this step also cor- 
rects for the time variation and vertex dependence of the 
detector response. The same method is used to deter- 
mine Evis for the MC sample. 

This analysis uses data collected by the Double Chooz 
far detector between April 13, 2011 and March 15, 2012, 
which is the same time-period used in the latest Dou- 
ble Chooz Gd analysis [6[. The total live time is 240.1 
days, which is different from 227.9 days used in the Gd 
analysis Q because of different analysis cuts. 

The IBD candidate selection is performed via the fol- 
lowing procedure. To reduce muon-induced backgrounds, 
we reject all events that occur less than 1 ms after a 
cosmic muon crosses the IV or the ID. We use PMT 
charge isotropy and PMT pulse simultaneity cuts to re- 
duce backgrounds caused by light emitted from PMT 
bases ("light noise") [1]. We apply the following coin- 
cident selection cuts to the remaining events: 0.7 McV < 
Eprompt < 12.2 MeV; 1.5 MeV < Edoiay < 3.0 MeV; 
10 < At < 600 /is; Ar < 90 cm. Furthermore, we 
reject prompt candidates that are coincident with a sig- 
nal detected in the OV. This veto, along with the 10 fis 
lower bound of the At cut, renders backgrounds due to 
stopped muons negligible. Finally, we apply a multiplic- 
ity cut to reduce fast neutron backgrounds. This cut de- 
mands that no trigger occur in the 600 ^s preceding the 
prompt candidate and that no trigger other than the de- 
layed candidate occur in the 1000 fis following the prompt 



candidate. 

The selection cuts yield a total of 36284 events. Among 
these IBD candidates are backgrounds due to uncorre- 
lated accidental coincidences, fast neutrons produced by 
muons traversing the nearby rock, long-lived cosmogcnic 
isotopes (mainly ^Li), and a small contribution from light 
noise. Accidentals are the dominant background, com- 
prising almost half the IBD candidate sample. 

We measure the rate and energy spectrum of acciden- 
tals by analyzing a sample of off-time coincidences. We 
collect this sample by looking for a delayed trigger be- 
tween 1 s -I- 10 ^s and 1 s -I- 600 ^s after a prompt candi- 
date event and applying a multiplicity cut for a period of 
1 s — 600 fxs to 1 s-l-1000 fis. To increase sample statistics, 
we open 124 consecutive windows after this first window, 
thus sampling accidentals between 1 s and 1.2 s after 
each prompt candidate. After correcting for inefficien- 
cies associated with this selection method, we obtain an 
accidentals rate of 73.45 ± 0.16 events/day. The result is 
cross-checked among multiple independent methods, and 
the quoted value includes the largest systematic uncer- 
tainty among them. 

The fast neutron background consists of a proton re- 
coil, the prompt event, in coincidence with the capture 
of the neutron, the delayed event. A single muon passing 
close to the detector may generate one or more fast neu- 
trons which traverse the IV and ID. We tag the number 
of IBD candidates in which fast neutrons are recorded 
simultaneously in the IV and ID by requiring > 2 IV 
PMT hits and an ID-IV pulse-timing correlation. We es- 
timate the tagging efficiency from an event sample with 
Eprompt > 12 MeV, following the same method as used 
for the Gd analysis Q. From this sample we obtain a 
spectrum shape and, using the tagging efficiency and 
sample purity, we calculate the fast neutron rate to be 
2.50 ± 0.47 events/day. 

Muon-induced radioactive isotopes which emit a neu- 
tron immediately following /3-decay, such as ^Li, can be 
a background to IBD reactions. As the lifetime of ^Li is 
257 ms, we use the correlation of the ^Li decay events 
to previously detected muons to estimate the ^Li back- 
ground rate. To increase the purity of ^Li in our sample, 
we consider only the subset of IBD candidates for which 
the spatial separation between the prompt event and the 
reconstructed muon track is within a defined distance. 
While ID PMTs are used to reconstruct the muon tracks 
in the Gd analysis 0, IV PMTs are used in this analy- 
sis to account for muons going through non-scintillating 
buffer liquid. To estimate the ^Li rate in this subsample, 
we fit the time difference At^ between the IBD candidate 
prompt events and preceding muons with an exponen- 
tial function characterized by the ^Li lifetime, plus a flat 
function to accommodate remaining accidentals and IBD 
candidates. The estimated rate is found to be consistent 
with that in the Gd analysis Q, accounting for the dif- 
ferent flducial volumes and selection efficiencies, and the 
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TABLE II. Summary of signal and background normalization 
uncertainties relative to the predicted signal. 
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FIG. 1. (Color online) Stacked histogram showing the prompt 
energy spectrum of neutrino candidates without background 
subtraction (black data points with statistical error bars). 
The red (grey) line is the best fit oscillation hypothesis. 
Also shown are contributions from accidentals (blue cross- 
hatched), ®Li (green vertical lines), fast neutrons (purple di- 
agonal lines), and correlated light noise (orange horizontal 
lines) . 



TABLE I. Summary of the number of observed IBD candi- 
dates and the predictions for the signal and background con- 
tributions used as input for the oscillation fit analysis. 



Source 


Predicted/Observed events 


h's prediction (no osc.) 


17690 


Accidentals 


17630 


Cosmogenic isotopes 


680 


Fast neutrons 


600 


Light noise 


80 


Total prediction 


36680 


Observed IBD candidates 


36284 



difference is included in the systematic uncertainty. We 
find a ^Li rate of 2.8 ± 1.2 events/day. Muon track re- 
construction efficiency is evaluated by a MC study and 
added into the systematic uncertainty. We estimate the 
shape and associated systematic uncertainty from MC, 
as was done in the Gd analysis 0. 

Finally, we found a small number of light noise events 
creating two consecutive triggers that are identified as 
IBD candidates. A volume cut on the reconstructed ver- 
tex is used to quantify the rate and Epiompt spectrum 
shape for this type of background. We estimate this back- 
ground rate as 0.32 ± 0.07 events/day. 

In this analysis, neutron detection efficiency e,! includes 
both the efficiency of the IBD selection and the fraction 
of neutron captures which occur on hydrogen. We eval- 
uate En from ^^^Cf neutron source calibration data and 
find it to be = 0.0846 ± 0.0018 in the target and 



Source 



Uncertainty 



Reactor flux 
Statistics 

Accidental background 

Cosmogenic isotope background 

Fast neutrons 

Light noise 

Energy scale 

Efficiency 



1.8 
1.1 
0.2 
1.6 
0.6 
0.1 
0.3 
1.6 



Total 



3.1 



e„ = 0.7853 ± 0.0036 in the 7-catcher. Weighting by the 
fraction of predicted IBD candidates in each region, we 
estimate the uncertainty in the detection efficiency over 
the entire fiducial volume as 1.0 %. Finally, we find that 
an uncertainty of 1.2 % accounts for the MC modeling of 
neutron migration, called spill-in/out between detec- 
tor subvolumes. Adding these factors in quadrature, we 
obtain a total detection efficiency uncertainty of 1.6 %. 

Energy scale uncertainty arises from three sources: 
time variation, non-linearity, and non-uniformity in the 
detector response. We treat the first two effects exactly 
as in Gd analysis Q . The third effect has a larger impact 
on the hydrogen analysis because of its extended fiducial 
volume. We estimate it by comparing data and MC from 
calibration source deployments in the 7-catcher. In total, 
we find an energy scale uncertainty of 1.7 %, as compared 
to 1.1 % used in the Gd analysis Q. 

The reference Epiompt spectrum is selected from the 
same reactor power-based MC sample generated for 
the Gd analysis 0. Systematic uncertainties on the ref- 
erence spectrum are the same as for the Gd analysis. We 
use the Bugey4 measurement to minimize the systematic 
uncertainty on the reactor neutrino flux prediction (gI-IT^. 
which is the dominant uncertainty in this analysis. The 
no-oscillation expectation for the number of neutrino can- 
didates is 36680 ± 520, including background. The pre- 
dicted number of events for both signal and backgrounds 
arc summarized in Tab. |T1 and uncertainties relative to 
the predicted signal statistics are shown in Tab. [Ill 

To extract sin^20i3 we compare both the rate and 
shape of the data to the reference Eprompt spectrum in 
31 variably sized energy bins from 0.7 to 12.2 McV. 
The fit procedure is identical to that used in the Gd 
analysis 0, Q, except that we use a single integration 
period and include the Ar cut efficiency as an addi- 
tional source of uncertainty. We use the MINOS value of 
Am2 = (2.32 ± 0.12) x lO^^eV^ as input for the fit 0. 
We find a best fit of 

sin226'i3 = 0.097 ± 0.034 (stat.) ± 0.034 (syst.) 

with xVDOF of 38.9/30. As in the Gd analysis d, 
we define statistical error as the portion of the 1 a 
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FIG. 2. (Color online) Top: Background-subtracted data 
(black points with statistical error bars) are superimposed on 
the prompt energy spectra expected in the case of no oscilla- 
tions (dashed blue line) and for our best fit sin^26i3 (solid red 
line). The best fit has xVDOF of 38.9/30. Solid gold bands 
indicate systematic errors in each bin. Middle: The ratio of 
data to the no-oscillation prediction (black points with sta- 
tistical error bars) is superimposed on the expected ratio in 
the case of no oscillations (blue dashed line) and for our best 
fit sin^2Si3 (solid red line). Gold bands indicate systematic 
errors in each bin. Bottom: The difference between data and 
the no-oscillation prediction is shown in the same style as the 
ratio (above). 



error which can be improved by collecting more data. 
This includes uncertainty from our current statistics (see 
Tab. [TT| and uncertainty on background shapes. We de- 
fine systematic error as the uncertainty which cannot be 
reduced simply by collecting more data. Figure [1] shows 
the complete spectrum of IBD candidates with the fit- 
ted background contributions, while Fig. [2] shows the 
background-subtracted Eprompt spectrum along with the 
best fit. The pull parameters from the fit are summa- 
rized in Tab. IIIII together with the input values. We 
have performed a frequentist study to determine the 
compatibility of the data and the no-oscillation hypoth- 
esis. Based on a Ax^ statistic, defined as the differ- 
ence between the at the best fit and at sin^20i3 = 0, 
the data exclude the no-oscillation hypothesis at 97.4 % 
(2.0 cr). A fit incorporating only the rate information 
yields sin226'i3 = 0.044 ± 0.022 (stat.) ± 0.056 (syst.). 
A simple ratio of observed to expected signal statistics 
yields R = 0.978 ± 0.011 (stat.) ± 0.029 (syst.) at the far 
site. 



TABLE III. Summary of pull parameters in the oscillation fit. 
The input values are determined by measurements, and the 
best-fit values are outcome of oscillation fit. 



Pull parameter 


Initial value 


Best-fit value 


Cosmogenic isotope [day~^] 


2.8 ± 1.2 


3.9 ±0.6 


Fast neutrons [day^^] 


2.5 ±0.5 


2.6 ±0.4 


Energy scale 


1.00 ±0.02 


0.99 ±0.01 


Am2(10"='eV^) 


2.32 ±0.12 


2.31 ±0.12 



The smaller best- fit value of sin^20i3 by the rate-only 
analysis can be explained by the ^Li background. The fit 
to the energy spectrum indicates a larger ^Li background 
contamination than the original estimate, although it is 
consistent within the systematic uncertainty. 

In summary, due to the low level of backgrounds 
achieved in the Double Chooz detector, we have made the 
first measurement of sin^20i3 using the capture of IBD 
neutrons on hydrogen. This technique enabled us to use 
a different data set with partially different systematic un- 
certainties than that used in the standard Gd analysis . 
An analysis based on rate and spectral shape information 
yields sin^26'i3 = 0.097 ± 0.034 (stat.) ± 0.034 (syst.), 
which is in good agreement with the result of the Gd anal- 
ysis sin226'i3 = 0.109 ± 0.030 (stat.) ± 0.025 (syst.) ji]. 
With increased statistics and a precise evaluation of the 
correlation of the systematic uncertainties, a combination 
of the two results is foreseen for the future. 
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